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FIG. 1: Two partly overlapping chains in a spherical cylinder 
of width D. The chains consist of TV beads of size a each. 
L denotes the chain extension, and R C 2c the distance of the 
centers of mass. 



Suckjoon Jun 

Faculte de Medecine, INSERM Site Necker, U571, 156 rue de Vaugirard, 75015 Paris, Franc$\ 

(Dated: February 1, 2008) 

We report molecular dynamics simulations of the segregation of two overlapping chains in cylindri- 
cal confinement. We find that the entropic repulsion between the chains can be sufficiently strong to 
cause segregation on a time scale that is short compared to the one for diffusion. This result implies 
that entropic driving forces are sufficiently strong to cause rapid bacterial chromosome segregation. 



I. INTRODUCTION 

Confined polymers play an important role in many 
industrial processes and biological systems. Examples 
range from membrane filtration and oil recovery to gel 
electrophoresis and protein translocation [USUI- Impor- 
tantly, recent technological development in nano- /micro- 
fluidics has made it possible to manipulate and trap 
biomolecules such as double-stranded (ds) DNA in con- 
fined environments with a characteristic lengthscale that 
is much smaller than the radius of gyration of the poly- 
mers [1, H, Also under biological conditions, DNA is 
often strongly confined, e.g. packed into a viral capsid Q, 
bacteria Q or the eukaryotic cell nucleus @ . 

In this article, we report molecular dynamics simu- 
lations that allow us to determine the typical speed of 
the segregation of initially mixed polymers in cylindrical 
confinement. This problem has particular relevance for 
the understanding of chromosome segregation in bacte- 
ria, where the nature of its underlying mechanism is cur- 
rently under debate. Here, the basic issue is whether the 
major driving force for segregation of duplicating chro- 
mosomes in strong confinement is physical (driven by en- 
tropy or mechanical "pushing" ) [1, [Trjl o r biological (such 
as cytoskeletal and motor proteins) [ill. Il2| . 

Our results show that the effective repulsion between 
two chains in a cylindrical geometry of confinement can 
be very strong. Typically, the segregation requires a time 
proportional to N 2 , which is much faster than the N 3 
timescale of chain diffusion, where N is the chain length. 
This suggests that for filamentous bacteria such as Strep- 
tomyces coelicolor [131 ] or cyanobacterium Anabaena [141 ] . 
the main driving force of chromosome segregation might 
be entropic and any additional mechanisms are for "op- 
timization." As we shall discuss later, our proposal is 
fully consistent with the recent results that chromosome 
segregation in some filamentous bacteria is a random pro- 
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II. THEORY 

Consider two linear chains with excluded- volume inter- 
actions, which are initially intermingled and confined in 
an infinitely long cylinder with a diameter D that is much 
smaller than the radius of gyration R g of the unconfined 
chains (Fig. [I}. As the two chains can gain conforma- 
tional entropy by demixing, they effectively repel each 
other. Note that the free energy cost for simultaneous 
overlap of n chains scales as n 9 / 4 in the dilute regime 
and then increases faster as n 3 in the more concentrated 
regime, independent of the chain length [l5|, [l6j]. Thus, 
for two intermingling chains illustrated in Fig. [TJ each 
blob-blob overlap contributes ~ fcsT (n = 2) to the 
free-energy cost. The potential of mean force for seg- 
regation between the two chains then is proportional to 
the total number of overlapping blobs, i.e., /3J r (i? C 2c) — 
(L oq — R C 2c)/D, where L eq is the equilibrium length of 
an isolated individual chain in the pore, and i? C 2c the 
center-to-center distance between the two chains. The 
effective repulsive force is then obtained as 
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and, thus, the equation of motion for the center of mass 
is 
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where M = Nm and T are the total mass and the ef- 
fective friction of the chain, respectively (to is the mass 
of a single monomer). Ignoring hydrodynamic interac- 
tions between monomers, one can assume that the fric- 
tions 7 on the individual monomers are additive, i.e., 
also r = 7V7. Then, the solution of Eq. [2] with an initial 
condition V C 2 C — at t — can be obtained as 



V c2c (t) 
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where r* = 771/7 is the "inertial" timescale. In practice, 
t ^> r* and hence the characteristic segregation speed is 
constant and given by 

k B T 
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The equilibrium length of confined chains, L oq , is pro- 
portional to N. Therefore, the time for reaching complete 
segregation, R c2c 



L eq , scales as 



t. 



L cq /V c ; 



N 
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This is time is much shorter than idiffi the typical time 
it takes a single chain to diffuse over a distance equal to 
its own length: 



^diff 
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However, the above considerations do apply for the ini- 
tial situation of complete overlap, R C 2 C ~ 0. In this case 
the system is in a state of unstable equilibrium, since the 
effective segregation force is F c g w 0. Hence, the sys- 
tem will initially show purely diffusive behaviour until a 
certain separation, typically R C 2 C ~ D, is reached. We 
refer to the time until segregation sets in as the "induc- 
tion time" U that should scale as N 3 . With increasing 
D, diffusion becomes easier, because the monomer con- 
centration decreases, and U decreases, while t s increases 
with D. Below, we show that for all practically relevant 
diameters, the segregation process is rate limiting. 

For real bacteria, entropic segregation already sets in 
during replication. Therefore, segregation always takes 
place for comparatively short pieces of DNA, so that the 
induction time does not play a role. 



III. SIMULATION METHOD 

In the molecular dynamics simulations, we model 
the polymers using a bead-spring model in a cylindri- 
cal compartment of diameter D; each chain consists of 
N beads of diameter a. The bead-bead and bead- 
compartment interactions were modeled by a Weeks- 
Chandler- Andersen potential (WCA) [jjj], which corre- 
sponds to the repulsive part of the Lennard- Jones poten- 
tial: 



UwCA{r) = £wca 
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TABLE I: The simulation parameters for the different runs. 
The first three tables give for different chain lengths N the 
simulated pore diameters D and the corresponding equilib- 
rium end — to — end distances L cq of a single chain. The 
last table contains the number T warm of timesteps used for 
equilibration of the interconnected chains, the number of 
timesteps T CO nfig between recorded configurations, and the 
number A con fl g of independent simulations runs with different 
random seeds. 



for r < \pla and elsewhere, r denotes the distance be- 
tween two bead centers for the bead-bead interactions, 
and the distance between the bead center and the com- 
partment minus a for the bead-compartment interac- 
tions. At r = a, the interaction energy is cwca = lfcsT; 
since the potential is quite steep, r will typically stay 
above 0.9er. This models soft beads of diameter a, whose 
centers cannot come much closer than a to each other, 
and cannot penetrate the wall (i. e. the wall imposes a 
constraint on the sphere centers, as depicted in Fig. []}. 
In the simulation, a defines the basic length scale and 
ewcA the energy scale. Our unit of mass is given by to, 
the mass of a bead. We choose the temperature such 
that ksT/e = 1. Having specified our basic units, the 
time unit is given by twca = cr a/to/ cwca = 1. In the 
following, we will omit these units. 

The springs between the beads in a chain were formed 
by the FENE (finite extensible nonlinear elastic) poten- 
tial 
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where r is the distance of the bead centers, rp is the 
radius at which the potential becomes singular, and ep 
is the interaction strength. In the present simulations, 
we chose ep = 10 and tf = 2. In combination with the 
WCA potential this results in a typical bond length of 
1.027. 

We simulate this system using the simulation package 
ESPResSo [18| . To propagate the system, we employ 
a velocity- Verlet MD integrator with a fixed time step 
of 0.01; the system is kept at constant temperature by 
means of a Langevin thermostat with a fixed friction of 
7 = tot^J-, a , so that t* — twca = !• Other parameters 
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FIG. 2: Example simulation run for D = 7, N — 200, starting 
from the removal of the interconnecting bonds. The two gray 
bands give the total extents along the tube axis of the two 
chains for a single run, the black lines the positions of their 
centers of masses. The positions are relative to the total cen- 
ter of mass R of the system and rescaled by the equilibrium 
length L eq of a single confined chain. 



vary for the different simulation runs, see Table Q] Our 
simulation procedure contains four steps: 

The system is initially prepared in a "ladder" configu- 
ration formed by two interconnected zig-zag strands, i.e., 
the system consists of two linear chains where the i— th 
bead of one chain is bonded to the i-th bead of the other 
chain, in addition to the bonds to its neighbors within 
the same chain. 

To equilibrate the system, we simulate for T warm steps 
with a "soft" WCA potential, i.e. a WCA potential that 
has been modified such that the potential is linear for 
distances smaller than a radius rf c . We reduce rf c gradu- 
ally during the equilibration phase, so that the potential 
converges to the plain WCA interaction. This proce- 
dure allows more overlap between beads during the ini- 
tial equilibration phase, which helps the bonds to quickly 
relax to the equilibrium length. 

After the equilibration of the interconnected chains, 
we remove the interconnecting bonds to obtain two sep- 
arate chains whose centers of mass very nearly coincide. 
The chains are stretched by about 10 — 20% compared 
to a single chain in confinement due to the cross linking; 
however, their length relaxes quickly to almost the same 
length as a single chain once the cross linking is released. 
The timescale for this relaxation is negligible compared 
to the segregation time. 

We continue to simulate the system until the two 
chains have segregated, i.e., until the chains do not over- 
lap and their centers of mass are separated by at least the 
equilibrium length L cq of a single chain, which had been 
determined beforehand by separate simulations. During 
this run, we record configurations every T con gg simulation 
steps. 
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FIG. 3: Schematic view of the segregation process. A line is 
fitted to the center of mass distance R C 2c in the range from D 
to L eq — D. From this fit, the induction time U is determined 
as intersection with R C 2c = 0, and the segregation speed V C 2c 
as its slope. The segregation time is then t s = L fxl /Vc2c- 

This procedure is repeated N coa Qg times (see Table HJ, 
resulting in iV con fig independent data sets similar to 
Fig. [21 For each of these data sets, we calculate the dis- 
tance Rc2c(t) of the centers of mass of the two chains 
parallel to the cylindrical compartment as a function of 
time. Initially, R C 2c is zero due to the preparation of 
the system, and stays close to zero during the induction 
time. Eventually, segregation sets in, and i? C 2c grows 
rapidly until i? C 2c = L cq is reached, at which time the 
chains do not overlap anymore. Further increase in R C 2 C 
is only due to diffusion and is therefore much slower. 

Fig. [3] displays schematically how we extract the in- 
duction and segregation times from each run: we fit a 
linear function (t — ti)V C 2 C to the range in which R C 2c(t) 
is between D and L cq — D. Here, U is the extrapo- 
lated onset time of segregation and V C 2 C is the speed 
with which the two centers of mass separate in the linear 
regime. We always find linear segregation behavior for 
D < i? C 2c < L c(i — D. The lower limit implies that the 
chains are separated by at least one blob diameter, the 
upper limit guarantees that there is at least one blob-size 
overlap left. 



IV. RESULTS AND DISCUSSION 

As can be seen in Fig. [4] our simulation clearly sup- 
port the scaling prediction V C 2 C ~ 1/(ND) (Eqn. Q). 
The prediction that the segregation time, scales as N 2 
is only recovered for small tube diameters. This is not 
unexpected, because when D approaches L eq , the simple 
blob prediction for L oq breaks down [l9[ and the segrega- 
tion time levels off at the relaxation of a free chain. The 
segregation speed relation Eqn. ((4]) however seems to be 
quite robust even for finite systems. 

The measured average induction time ti is shown in 
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FIG. 4: Measured segregation speed V C 2c rescaled by 1/N. For 
better visualization, we actually plot its inverse, l/(iW C 2c). 
The dashed line demonstrates the linear scaling of the segre- 
gation time with D. The inset shows the segregation time t s 
rescaled by iV 2 , as a function of D. 



FIG. 6: Monomer density p(r) along the cylinder-axis for 
N = 200, D — 8, averaged over all configurations with a 
center of mass distance R C 2c ~ (a), R C 2c = 20 (b), i? C 2c = 48 
(c) and R C 2c = 80 (d). The graphs are centered around the 
systems center of mass R. 
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FIG. 5: Measured induction times i, rescaled by TV 3 . The 
lines are just guides to the eye. 



effective entropic force F e g. The overwhelming part of 
the chain "demixing" is due to directed segregation (see 
Fig. In other words: it may take a while for the 
system to start segregating, but the segregation process 
itself is always governed by the effective entropic force. 

Fig. [S] shows the average monomer densities of the two 
polymers for different center-of-mass distances i? C 2c- As 
predicted, the monomer densities of each chain in the 
overlap region are almost unaffected by the presence of 
a second polymer in the same space. Hence, the initial 
monomer density is almost twice as large as for a sin- 
gle chain. During segregation, the monomer densities 
of the individual polymers increase somewhat. In fact, 
the snapshots show that the polymers have very nearly 
separated at i? C 2 C = 48, which is significantly less than 
L cq = 61.1; this demonstrates that the polymers deform 
during segregation: the entropic driving force is strong 
enough to compress the polymers. After demixing, the 
chains expand to their equilibrium length. 



Fig. [5] The distribution of the induction times has a 
long tail, which makes it difficult to sample t% accurately. 
Keeping this caveat in mind, we find that for small tube 
diameters the ti's computed for different chain lengths 
can be made to collapse if we assume A^ 3 -scaling, as ex- 
pected for a diffusive process. Moreover, we do observe 
the expected decrease of U with increasing D. For larger 
tube diameters, the induction time increases again; this 
is probably due to the fact that for larger D the segre- 
gation and induction times cannot be clearly separated 
(ti/t s = 0(1)) for larger diameters. In fact, in our sim- 
ulations, the induction time seems to converge to about 
one quarter of the segregation time for all N and D <J 4. 

We stress that for highly confined chains, the diffusive 
process is only responsible for the segregation over the 
tiny initial separation necessary to obtain a significant 



V. CONCLUSIONS 

Our simulations support the scaling prediction that the 
entropically driven segregation of two confined chains re- 
quires a time proportional to iV 2 . For long chains, this 
time is much shorter than the diffusive time that scales 
as iV 3 . We stress that this speed up of entropically 
driven segregation does not involve any active (energy- 
consuming) process. Considering the geometry of con- 
finement and the length scales of (filamenteous) bacteria, 
our results strongly suggest that the partitioning of dupli- 
cated chromosomes in these organisms is, at least partly, 
entropy-driven. Since the segregation sets already in dur- 
ing replication, there is no initial "induction" regime for 
bacteria. Indeed, the recent data obtained by Hu et al. 
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on the cyanobacterium Anabaena sp. PCC 7120 suggests 
that MreB, a bacterial actin homologue that is specu- 
lated by some as a "track" for transporting chromosome 
by putative motor proteins, is important for cell shape 
but not for chromosome segregation [14|]. Moreover, they 
also have shown that the ratios of DNA content in two 
daughter cells have a much wider distribution than in the 
case that the two cells were identical. This suggests that 
chromosome partitioning is a random process, in good 
agreement with our entropy-driven segregation process 
we presented in this article. 
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